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Abstract. We consider a stabilization problem for a piezoelectric system. We prove an expo- 
^-^' nential stability result under some Lions geometric condition. Our method is based on an identity 

C^ , with multipliers that allows to show an appropriate observability estimate. 
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1 Introduction 



< 

i -p, ' We consider the dynamical behavior of a piezoelectric system (which means the ability 

C^ , of some materials, like ceramics and quartz, to generate an electric field in response to 

applied mechanical stress), where a proper modeling involves the displacement vector, the 
electric field and the magnetic field, which are governed by the elasticity system coupled 
with Maxwell's equations. This system plays an important role in various applications in 
►^ , structural mechanics and in mechatronics, for such a model we refer to |10[ 115 ]. 

^*7^ . Let 17 be a bounded domain of IR'^ with a Lipschitz boundary T. In that domain we 

1^ ' consider the non-stationary piezoelectric system that consists in a coupling between the 

p^ . elasticity system with the Maxwell equation. More precisely we analyze the partial dif- 

• ' ferential equations based on the following relations between the stress tensor, the electric 

^— V . displacement and the magnetic induction: 



(1.1) a.ij{u,E) = aijkijkiiu) - CkijEk Vi, j = 1,2,3, 



k> ; (1.2) A = eij£^j+e»fcz7fei(M)Vi = l,2,3 



(1.3) B = fiH. 

The equations of equilibrium are 

(1.4) d^u, = djaj,Wi = 1,2,^ 
for the elastic displacement and 

(1.5) dtD = cm-lH, dtB = -curlS 

for the electric/magnetic fields. 

This system models the coupling between Maxwell's system and the elastic one, in which 
E{x,t), H(x,t) are the electric and magnetic fields at the point x € fl at time i, u(x,t) is 
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the displacement field at the point a: G 51 at time t, and 'yij{u)^ ^-^ is the strain tensor given 

by 

, . 1 f dui duj 

Here a = {o'ij)i j=i,D = (_Di,_D2, £'3), and B = {Bi,B2,B^) are the stress tensor, electric 
displacement, and magnetic induction, respectively, e, /i are the electric permittivity and 
magnetic permeability, respectively, and we will assume that they are positive real numbers. 
The elasticity tensor {aijki)i,j.k,i=i.2,3, is made of constant entries such that 

O-ijkl = ajikl = aklij 

and satisfies the ellipticity condition 

(1-6) aijkili]lki > o-olijlij, 

for every symmetric tensor (7ij) and some ao > 0. The piezoelectric tensor Ckij is also made 
of constant entries such that 

^kij ^kji • 

For shortness in the remainder of the paper introduce the tensor a{u) — {aijki^ki{u))1 ^^i 
and let Vct be the vector field defined by 

while for a tensor 7 = (7ij)^,-^i, and a vector F — (-Fi, -P^j -F3), we set 

67 = {.eiMlki)l=i e^F = (e^kiF^)l^i^^. 

These last notations mean that e corresponds to a linear mapping from M'^^^ 
that e^ is its adjoint. With these notations, we see that (|l.ip is equivalent to 

a{u,E) =(t{u) -e^E, 



into 



and 



while (|1.2[) is equivalent to 



D ^ eE + e-f{u) 



The system (|l.ip - (ll.3p is completed with the boundary and Cauchy conditions. This 
means that we are considering the following system 



(1.7) 



dtu 



2- Vcr{u,E) =0 in Q ■.= nx]0,+oo[ 
dtD - cm\H = in Q, 



fidtH + cuilE = in Q, 

div(L>) = div(/iiJ) = in g, 

H xu- {Q*dtu) xiy+{E X I') X 1^ = on j: = r X {0, +00), 

a{u, E) ■ u + Q{E xv)+ Au + dtu = on S, 

u{0) = uq, dtu{0) = ui in 17, 

E{0) =. Eo, H{0) = Ho in n, 

where v is the unit normal vector of dfl pointing towards the exterior of H., A is a. positive 
constant and Q is a function from F into the set of 3 x 3 matrices with the regularity 

QeL°°(F,C3x3). 



Remark 1.1 Note that the image of Q of normal vector fields plays no role in the boundary 
conditions appearing in (|f .7p . Indeed for X G C'^, let X^ — [X ■ v)v and X.,- = X — X^ be 
the normal and tangential components of X respectively, by writing 

QyX = QX^, QrX = QXt, 

we get the splitting 

But by definition QIX is orthogonal to the tangent plane. Therefore Q{E xv) ~ Qt{E x v) 
and {Q*dtu) x v = {Q*^dtu) x v, which means that the normal part Qi, of Q does not 
contribute to the boundary conditions. 

Boundary or internal stability of the second order elliptic systems, like the wave equation 
or the elasticity system, have been studied by many authors, let us quote O [51 [H [TH 
[T4l [21] [23l [24l [27l [28] among others. Similar results for Maxwell's system can be found 
in [T6l[20l[22l[2lllll[3Tl[25l[34]. The combination of these results to the piezoelectric 
system in some particular cases has been treated in [17l [TH [32]. For the quasi-static case 
(corresponding to the hypothesis that E is curl free, hence the gradient of a potentiel), we 
can refer to [Tgi[^ [^. 

In [T5] the authors consider the above problem in the case Q — Q with eventually dis- 
continuous coefficients and an additional memory term and prove the exponential decay 
rate of the energy if A is small enough and if fi satisfies some geometrical conditions (star 
like shape). On the contrary in [32] j the author treats the case Q = I and some nonlinear 
feedback terms, but with the choice of e such that V(e^£^) = £^cui\E for some real number 
^ (case excluding the natural condition Ckij = Skji) and proves the exponential decay rate of 
the energy in the case of linear feedbacks if H. is strictly star shaped with respect to a point. 
In that last paper the author combines the multiplier technique with the one from [4 , where 
the authors uses some tangential integration by parts and a technique from [6 . Our goal is 
here to perform the same analysis for the general system p.7p . For Q e L°°{r,C'^^^) we 
prove that the system is well-posed using semigroup theory. On the other hand using the 
multiplier method (see [T^) and a technique inspired from [21 [H [TSl [32] to absorb a zero 
order boundary term, we show that the system is exponential stable if Q = a / for some 
scalar continuously differentiable function a such that Va is small enough. 

The paper is organized as follows. The second section deals with the well-posedness of the 
problem. In the last section we give the main result of this paper which is the exponential 
stability of the piezoelectric system and its proof. 

2 Well-posedness of the problem 

We start this section with the well-posedness of problem (|1.7p . At the end we will check the 
dissipativeness of p. 71) . 

Let us introduce the Hilbert spaces (see e.g. |25[l30j ) 

(2.8) J{n) ^{E e L^{n)^\divE = in fl}, 

(2.9) n = H^inf X L^{nf x L^{fif x j{n), 

equipped with the norm induced by the inner product 

{E,E'), = I eE{x) ■E'{x)dx,\fE,E' G J{n), 

((w, V, E, H), iu\ v\ E\ H'))^ = (u, u')i + {v, v')o 

+ {E, E'), + {H, H')^y{u, V, E, H), (u', v' , E' , H') e %, 



where we have set 

{u,u')o— I u{x) ■ u' [x) dx , 



{u,u')i— I (t{u){x) : '-f{u'){x)dx + A u{x) ■ u'{x)dS, 
Jn Jr 

with the notation 

cr{v) :7(u') := cry(w)7y(u')- 

Now define the hnear operator A from H into itself as follows: 

(2.10) D{A) ^ {{u, V, E, H) £ H| Vf7(u, E), curl^, curlff e L^{nf;v e H\nf; 

E X i/,H X u e L^(r)^ satisfying 

(2.11) H xiy-{Q*v) xiy + E xiyxiy = Oonr, 

(2.12) cr{u, E)-iy + Au + v + Q{E x v) = Q on T}. 

For all (u, w, E, H) e D{A) we take 

A{u,v,E,H) ^ {v,V(T{u,E),e~^{cnr\H - e-f{v)),~^i~^cm\E) . 

The boundary conditions (|2.11l) and (|2.12p are meaningful since for {u,v,E,H) G D{A), 
from section 2 of [3] the property \7a{u,E) e L^{il)^ implies that a{u,E) ■ v belongs to 
7f-i/2(p-)3_ gjj^pg i-jjg properties m, w G R^(yt)^ imply that Au + v belongs to iJi/2(r)^, the 
boundary condition (|2.12p has a meaning (in 7J~^/^(r)^) and furthermore yields cr(w, E)-v ^ 
L^(r)^ (because Q(E x v) ^ L^(r)^). Similarly the properties of H and v give a meaning 
to the boundary condition (|2.11l) (as an equality in I/^(r)'^). In summary both boundary 
conditions (|2.1ip and (|2.12p have to be understood as an equality in L^(r)'^. 

We now see that formally problem (I1.7P is equivalent to 

r ^- All 

when JJ = {u, dtu, E, H) and Uq = (uo, ui, Eq, Hq). 

We shall prove that this problem (|2.13p has a unique solution using semigroup theory by 
showing that ^ is a maximal dissipative operator. 

Lemma 2.1 A is a maximal dissipative operator. 

Proof: We start with the dissipativeness: 

iAU,U)n <0,yu eD{A). 
From the definition of A and the inner product in Ti, we have 

(AU, U)n = {v, u)i + (Va(u, E), v)o 

+ {E- (curliJ - e"f{v)) - curlE' • H} dx, 
Jn 

for any {u,v,E,H) G D{A). Lemma 2.2 of [31j and Green's formula yield equivalently 

{AU, U)h = («, u)i - f aiu, E) : j{v) dx 
Jn 

ej{v) ■ E dx 
n 

+ I {{a{u,E)-u)-v + {E X v)-H)dS, 



for any (m, v, E, H) G D{A). Using the definition of the inner product (•, •)! and the bound- 
ary conditions !|2.1ip and (|2.12p . we arrive at 



{AU, U)n^- {|f r + \Ex iy\^}dS < 0, 

for any {u,v,E,H) £ D{A). 

Let us now pass to the maximahty. This means that for at least one non negative real 
number A, A / — ^ has to be surjective. Let us show that indeed I — A is surjective. This 
means that for all (/, g, F, G) in Ti, we are looking for {u, v, E, H) in D{A) such that 

(2.14) {I-A){u.v,E,H)^{f,g,F,G). 

From the definition of A^ this equivalently means 

u-v = f, 
(o^K^ ) v-Va{u,E)= g, 

^ ^^' ^ E - e-^cmlH - ejiv)) = F, 

H + fi-^curlE = G. 

The first and fourth equations allow to eliminate H and v, since they are respectively 
equivalent to 

(2.16) v = u-f, 

(2.17) H = G- fi^^cuTlE. 

Substituting these expressions in the second and third equations yields formally 

(2.18) u-V(t{u,E) = f + g, 

(2.19) eE + curl(^"^curl£;) + e-/{u) ^ eF + curlG + e-f{f). 

This system in (m, E) will be uniquely defined by adding boundary conditions on u and 
E. Indeed using the identities p.l6p and (|2.17p . we see that (|2.1ip and (|2.12p are formally 
equivalent to 

(2.20) -^"^curli; X J/ + Q*u xv+{Exv)xv^-Gxv + Q* f x i' on T, 

(2.21) a{u, E)-i' + Au + u + Q{E x i^) = / on P. 

By formal integration by parts we remark that the variational formulation of the system 
([2T8| - ([2T9l) with the boundary conditions ([2:20| - ([23T|) is the following one: Find {u, E) eV 
such that 

(2.22) a{{u, E), {u' , E')) = F{u' , E'),^{u\ E') G V, 

where the Hilbert space V is given by V^ = H^{Q)^ x W when W is defined by 
W ^{E (^ L^{VLf\cm:\E e L^iflf and E xv e L^{Tf}, 

with the norm 

I \E\ 1^ = f i\E\^ + \cur\E\^)dx + [ \E x iy\^ dS, 
Jn Jr 

the form a is defined by 



a{{u,E),{u',E')) = {(j(u,E) ■.j{u')+u-u'}dx 

Jn 

+ / {^"^curlS • curl£;' + eE ■ E' + e-f{u) ■ E'} dx 



+ I {{Exv)- {E' xv) + {A+ l)u ■ u' + Q{E x v) ■ u' - Q*u ■ {E' x v)} dS, 



and finally the form F is defined by 

F{u\E')^ f {{f + g)-u' + {eF + ej{f))-E' + G-cmlE'}dx+ j {f -u' ~{Q* f xv)- E') dS. 
Jo, Jv 



10. Jv 

We easily see that the bilinear form a is coercive on V since 

a{{u,E),{u,E)) = / {<t(m) ; 7(m) + |Mp} dx 
Jo 



( {li-^\cnT\E\^ + e\E\^} dx 
Jo. 

f{\Exi^\^ + {A + l)\u\^}dS, 



which is clearly greater than ||w||^i(s^-)3 + ||-E'||^ by the ellipticity assumption on the elasticity 
tensor. Hence by the Lax-Milgram lemma, problem (|2.22p has a unique solution (u, E) e V. 

To end our proof we need to show that the solution (u, E) € V oi (|2.22p and v, H given 
respectively by (PTTB)) . (PT7| are such that (w, u, E, H) belongs to D{A) and satisfies (pji)) 
(or equivalently (|2.15p ). First taking test functions u' in 2?(51)'^ and E' = 0, we get 

Va{u,E)+v^g mV'{n). 

This implies the second identity in (12.151) as well as the regularity \7(t{u,E) G L^{il)^ (from 
the fact that v, curl_E as well as g belongs to that space). 

Second we take test functions u' ~ and E' — x with x £ 2?(il)^ by Lemma 2.3 of |31) 
we get 

eE ~ curliJ + e-^iu) = eF in Vi^l). 

This means that the third identity in (|2.15p holds as well as the regularity curli? £ _L^(J7)^. 

Thirdly taking test functions v' G H^i^Y ^^'^ E' = x with x G C°°(fj)'^ and applying 
Green's formula (see section 2 of [3] and Lemma 2.2 of f3T]), we get 



(cr(u, E) ■ u, v') - {H xiy)- E'dS+ {Q{E x v) ■ u' ~ {Q*u x v) ■ E' dS 
+ \ {{Exv)- [E' xv) + {A + l)u ■ u'} dS = 0. 



Jr 

This leads to the boundary conditions (|2.1ip and (|2.12p since u' (resp. x) was arbitrary in 
H^{^)^ (resp. in C°°{0,)^) whose trace belongs to a dense subspace of L^{r)^. 

Finally from (|2.17p and the fact that iiG is divergence free, fiH is also divergence free. 

■ 

Semigroup theory |331 [36] allows to conclude the following existence results: 

Corollary 2.2 For all {uq,ui, Eq, Hq) G H, the problem (J. 71) admits a unique (weak) so- 
lution {u,E,H) satisfying {u,dtu, E, H) G C(R+,'H), or equivalently u G C^{W+,L'^{Vl)'^) fl 
C(]R+, if 1(17)3), E G C(R+,L2(r2)3) and H G C(E+,J(fi)). If moreover {ua,ui, Eo,Ho) 
belongs to D{A) and satisfies 

div(e7(uo) + sEq) — in fl, 

then the problem ( tj.?)) admits a unique (strong) solution (u, _E, H) satisfying (m, dtu, E, H) G 
C^{R+,'H) r]C{R+,D{A)), or equivalently satisfying u G C^{R+, L^{n)^)nC^{R+, H^{nf), 
E eC^{R+,Jin))r]C{R+,W), H eC^{R+,Jin))r]C{R+,W), satisfying ^KW-^MIIM) for 
a.e. t (with v — dtu), as well as 

Va{u,E) GC(M+,i2(f7)3). 
6 



Note that, in that last case, D = e^{u) + eE satisfies in particular 

divD = m rj X M+. 

We finish this section by showing the dissipativeness of our system. 
Lemma 2.3 The energy 

(2.23) £{t) = I f {\dtu{x,t)\'' + c7{u){x,t):-f{u){x,t))dx + ^ [ \u{x,t)\^dSix) 

+ ^ f ie\£{x,t)\^+^i\Hix,t)\^)dx 
is non-increasing. Moreover for {uq, ui, Eq, Hq) G D(A), we have for all < S < T < oo 

(2.24) £{S)~S{T)^ f f{\E{x,t)xi^\^ + \dtu{x,t)\^}dSdt, 

Js Jr 

and for all t > 

(2.25) dt£{t) = - f{\E{x,t) X z/p + \dtu{x,t)\^} dS. 

Proof: Since D{A) is dense in H it suffices to show (|2.25p . For (uo,ui,Eo,Ho) G D{A), 
from the regularity of u, E, H, we have 

dt£{t) = f{dfu-dtu + a{u):j{dtu)}dx + Afdtu-udS 
Jn Jr 



{sE ■ dtE + nH ■ dtH} dx 

By 1131), we get 

dt£{t) = I {dtu-V(7{u,E) + a{u):j{dtu)}dx + A I dtu-udS 
Jn Jr 

+ f {E- {curlH - E-i{dtu)} dx 
Jn 

= [A {u{t),dtu{t),E{t), H{t)) , {u{t),dtu{t),E{t), H{t)))^ 



We conclude by Lemma [2. II ■ 

3 Exponential stability 

Lr this section we prove the main result of this paper, namely the exponential stability of our 
system ()1.7p when il is strictly star-shaped with respect to a point x^. This result is based 
on an identity with multipliers proved in [18] that allows to show the next observability 
estimate. 

Theorem 3.1 Assume that there exists xq G K" and (5 > such that 
(3.26) m{x) ■ v{x) >5 Vx G dVi, 



where m{x) — x — xq- Assume also that Q = al with a continuously differentiable function 
a from F to C. Set Ca = max|VQ:|. Let {u,E,H) be the strong solution of problem jj.?)) . 

Then there exists a positive constants C (independent of a) such that for all T > 0, and all 



9 , there exists a constant C{9) (independent of T) such that the next observability estimate 
holds: 

(3.27) T£{T)<{C{e){l + CaT) + eT)£{0) + C f {\dtu\^ + \Er\^) dSdt, 
where Et = T x (0,r). 

Proof: First the identity (3.9) of [TB] with to = and 1^9(2;) — \x — .to|^/2 yields 

(3.28) T£{T)^r+ I V{x,t) dS{x)dt, 

where we have set 

r = ~2 Idtu- {u + {m-V)u)} + fi{mx H)- {eE + e^{u)}\dx , 

V = 2{tdtu+{m-V)u + u}-cr{u,E)i^ + m-i^{\dtu\^ ~a{u) ■.j{u)+e\E\^ +H\H\'^} 
+ 2t{H y. E)-v ~ 2eE -vE -m- 2fiH -vH -m 
— 2(m X e^{u)) ■ {E x v). 

Using the boundary conditions from (|1.7p . we see that 

V = -2tdtu{Q{E xiy)+Au + dtu) + A 

+ m-iy{e\E\^ +ii\H\^} 

- 2t{Q*dtu) xv)- E- 2t\Er\'^ - 2eE^{E„m ■ v + Er ■ rur) - 2fiH^{H^m -v + Hr-mr 

— 2{m X e^{u)) ■ {E x v), 

where we recall that E^, — E ■ i/, Er = E — E^v and 

A = 2{(m • V)m + u] ■ a{u, E)v + m ■ v{\dtu\'^ - a{u) : 7(w)}. 

By Young's inequality, there exists C > such that for all /3i , /32 > 

V < ~2Atudtu 

- [m ■ u ~ (i2){e\E,\^ + ^l\H,\'') 

+ A + C{l + ^ + ^)\Exv\^ + C(l + U\H^ v\^ + A7(«) : ^[u). 

P2 Pi P2 

By using again the first boundary condition from p.7p . we get for all /3i, /32 > 

(3.29) V < -2Atudtu + Cil + ^)\dtu\^ + C{l + ^ + ^)\Exiy\^ 

P2 P2 Pi 

- (m • z. - /32)(e|^.P + Ml^.n + A + /3i7(u) : 7(")- 
Let us transform the first term of this right-hand side: 

-2A tudtudSdt= ~A t—u^dSdt, 
and by an integration by parts in time, we get 

T 




-2A / tudtu dSdt = A u^ dSdt ~ A I tu^ dx 



This proves that 

(3.30) -2A I tudtudSdt <A u^ dSdt. 



Let us now estimate the term A. First using the second boundary condition from (|1.7p . we 
see that 

A = -2{(to ■ V)u + u} ■ {Q{E xv)+ Au + dtu) + m ■ i'{\dtu\'^ - a{u) : 7(u)}. 

Using the elhpticity assumption (|1.6I) and condition (|3.26p we obtain 

(3.31) A < -2{{7n-\7)u + u}-{Q{Exiy)+Au + dtu) + m-i^\dtu\'^~aod-f{u):-/{u) 
< -2m ■ Q{E Y.V)- 2A\u\^ - 2u ■ dtu - 2(m ■ V)u ■ Q{E x i^) 
— 2A{m ■ V)ii • u — 2(to • \7)u ■ dtu + m ■ v\dtu\ — aQd^{u) : j{u). 

We need to estimate some terms of this right-hand side. First as before an integration by 
parts in time yields 



-2/ udtudSdt<A \u{x,t ^ 0)\^ dS{x) <2£{0). 

JT.T Jt 

As in [U [13], one can show that 

(3.32) -2yl/ {m-\7)u-udSdt < t- / \uf dSdt + Oi [ j{u) : j{u) dSdt, 

as well as 
(3.33) 

/ {m-V)u-dtudSdt<C£{0)+^ [ i\u\'^+\dtu\^)dSdt+e2 f -/(u) : -f{u)dSdt,\^ei,92 > 0. 

By Young's inequality we clearly have 

(3.34) / u-Q{Exv)dSdt<C I {\u\'^ + \Er\'^)dSdt. 

Now we notice that 

(m • V)m • Q{E xv)^ {Q*{m ■ V)u) ■ [E x z/), 

and for any k = 1,2,3, we may write 

{Q*{m-V)u)k = Qlj-rriidiUj 

= 2QljTnijij{u) - QljiTiidjUi 

= 2QljTnijij (u) + QljUidjirii - Qljdj (niiUi). 

The two first terms of this right-hand side will be estimated by Young's inequality and it 
therefore remains to estimate the last term, namely by the previous identities we have 



(3.35) / {m-\7)u-Q{Exi^)dSdt<e3 j{u) : -f{u) dSdt 

-Fc/ {\u\^ + {\ + ^)\Er\^)dSdt~ \ {Q*V{m-u))-iExiy)dSdt,ye3>0. 



'E 

Now using Green's formula, we see that 



/ {Q*W{m ■ u)) -{E xv) dSdt = f {curl(Q* V(m ■ u)) ■ E - Q*V{m ■ u) ■ cm\E} dSdt, 

where Qt — ^ x (0,r). Now using the fact that Q{x) — a{x)I, and that curlE' = ndfH, 
we obtain 

/ (Q* V(m • u)) -{Exv) dSdt = / {(Va x V(m ■ u)) ■ E - Q*V{m ■ u) ■ ^dtH} dSdt. 

J YIt ^ Qt 



For this last term, we first integrate by parts in time and get 
f Q*V{m-u)- ndtHdSdt = - f Q*V{m ■ dtu) ■ fiH dSdt+ f Q*V{m ■ u) ■ fiH dx . 

Jqt Jqt J^ ° 

An integration by parts in space leads to 
/ Q* V(m • u) ■ fidtH dSdt ^ / {Q*m ■ dtudw{nH) + m ■ dtuVa ■ (fiH)) dSdt 

Jqt Jqt 

- / Q*m-dtu{fiH)-iydSdt+ Q*\7{m-u)- fiHdx . 
These two identities and reminding that div(fj,H) = lead to 

{Q*V(m-u))-(EY.v)dSdt = / {{VoLY.V(rn ■ u)) ■ E - m ■ dtuVa ■ {[iB)] dSdt 

+ Q*m- dtuifiH) ■ V dSdt - / Q*\7{m ■ u) ■ fiH dx 



IT. 

By Young's inequality we find that 







/ {Q*V{m-u))-{Exv)dSdt<C{l + CccT)£{Q)+ [ {^\dtu\^ + 9i\H^\^} dSdt,yei > 0. 
This last estimate in (|3.35p leads to 



(3.36) / (m • V)u • Q{E x v) dSdt < C(l + c„T)£:(0) + 6j, I 7(u) : -f{u) dSdt 

+c f {|up + hdtuf + (1 + ^)\Er\^) + e4H,\^}dSdt,y03,ei > o. 

Now using again Young's inequality and the estimates (|3.32p , (|3.33[) , (I3.34p and (|3.36p into 
the identity p.3ip . we obtain that 



AdSdt < C{l + CaT)£{O) + {~ao6 + 9i+02 + 93) j{u) : "fiu)dSdt 

+c [ {{i + ^ + h\dtu\'' + {i + ^)\Er\^) + e4\H,f}dSdt 
+ c{{i + ^ + ^}f \u\^dSdtyei,92,e3,94>o, 

This estimate in p.29p and using p.30p . we get finally 

f VdSdt < C{l + c^T)£{0) + {-ao5 + /3i+9i + 92 + 93) f -f{u) : j{u) dSdt 

+ C f {I + ^ + I + -^ + ^)\dtu\' dSdt 

+ C f {1 + ^ + ^ + l)\Er\' dSdt 
JT.T Pi ^2 6*3 

+ C ( {{-m-v + l32)e\E,\'' + {{-m-u + p2)tJ^ + 9A)\H,\'']dSdt 
+ C{l + -3- + l}/ \u\^dSdt,y/3i,l32,ei,92,d3,9A>0. 

By choosing /3i , /?2 , ^i , ^2 , ^3 and 64 small enough, we have found that 
(3.37) / VdSdt<C{l + CaT)£{0) + C f i\u\^ + \dtu\^ + \Er\^)dSdt. 
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Coming back to (|3.28p and using again Young's and Korn's inequalities to estimate r, we 
obtain 

(3.38) T£{T)<C{l + CaT)£{0) + C [ {\u\^ + \dtu\^ + \Er\^)dSdt. 
Now invoking Lemma 13.71 below, we arrive at 

Te{T) < C(6')(l+c„T)£(0) + C / {\dtuf + \Erf)dSdt + e[ S{t)dt 

< {c{e){i + co,T) + eT)£{o) + c f {\dtuf + \Erf)dSdt,ye>o, 

reminding that the energy is non increasing. This is the requested estimate (|3.27p . ■ 

Remark 3.2 Note that the last term of the estimate p.35p is zero if Q* = Ql, but according 
to Remark II. 11 this assumption is meaningless. 

Theorem 3.3 Under the assumptions of the previous theorem and if c^ is small enough, 
there exist two positive constants M and lo such that 

(3.39) £{t) < Me~'^*£{0), 
for all strong solution {u,E,H) of il.T}) . 

Remark 3.4 The same method yields the same exponential stability result in the case 
where e, fi are positive functions satisfying some regularity and technical conditions. 

Proof: The estimate (|3.27p and Lemma [2731 vield 

T£{T) < (C(0)(1 + c„r) + eT)£{0) + C{£iO) - £(T)), V0 > 0, 
which may be equivalently written 

Now we choose 9 — ^ and Cq < J,j_, , with this choice p^y tends to C(^)cq + 5 < 

£{T) < r£{0). 



^ as T goes to infinity. Therefore for T large enough, we have found r G (0,1) such that 



Since our system is invariant by translation, standard arguments about uniform stabilization 
of hyperbolic system (see for instance [351131] ) yield the conclusion. ■ 

The key point in the above proof is to estimate appropriately the term L, [ul"^ dSdt in 
(|3.38p . Indeed a rough idea is to use the definition (|2.23p of the energy to get 

T 



C \u\^ dSdt < — / £{t) dt < — i-^(O) 
jet ^ Jo ^ 



Hence from the previous proof we obtain an exponential stability result only for A small 
enough (depending on a constant C that is not known explicitly, see nevertheless [2J). In 
order to prove the stability result for any positive A, we then need to estimate L, |u|^ dSdt 
in a different way. Its proof is based on the use of a solution z of a stationary problem (see 
[HI H [131 [32] and below) such that z = u on F. Muhiplying the first identity of ([TT7|) by z, 
integrating by parts and using the second boundary condition in (jl.7p . the term L |w|^ dSdt 
naturally appears. For standard problems (see JSJ [U [T31 [32] ) this term is estimated using 
elliptic regularity results on z. Here the specificity of our piezoelectric system requires a 
more careful analysis. We start with the stationary problem mentioned before. 
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Lemma 3.5 Let (m, E, H) be a strong solution of il.T]) . Then there exists (z, x) £ H^{il)'^ x 
-ffg(f2) (depending on t) weak solution of 

V(cr(z) - pJVx) = m r^, 

(3.40) { div(eVx + e7(z)) = Q mVt, 

z — u, X = onV . 

Moreover there exists a positive constant C (independent of t) such that 

(3.41) f \z\^dx<C f \u\^dS <^£{t), 

Jn Jr ^ 

(3.42) / \dtz\^ dx<C f \dtu\^ dS < -Cdt£{t). 
Jn Jt 

Proof: Inspired from [SJ [21 [131 [32] for each i > we consider the weak solution (z,x) 
(depending on t) of p.40p . This solution is characterized by z = w + w where (w, x) G V := 
HQ{il)^ X Hq{U) is the unique solution of 

(3.43) diiw, x), {w\x')) = -5((", 0), {w', x')Mw',x') e 1>, 
where 

Hiw,x),iw',x')) = / {('^(li') -e^Vx) : "/{w') + (eVx + e7(u;)) ■Vx}dx,y{w\x) € K 

The above problem has a unique solution since the bilinear form a is coercive on V (conse- 
quence of Korn's inequality). 

A direct consequence of p. 431) is that 

diiz,x)Aw\x')) = OMw',x')^V. 
By taking as test function w' = w — z ~ u and x' — X-, we find that 

a{{z,x),{z^x)) = a((z,x),(u,0)), 
which implies 

(3.44) / {a(z) : jiu) - e^ Vx : 7(")} dx = d{{z, x), (z, x)) > 0. 
Jn 

Note further that the coerciveness of a leads to 

lklli,a + llxlli,a<Ch||i,n, 
and then to 

(3.45) \\zh.n + \\x\\i.n<C\\uh,n<C£{ty^^ 

where ||u||s^o = ||u||i/=(o). 

Now we consider the adjoint problem; Find {'w*,x*) £ ^ solution of 

V((T(u;*) + e^Vx*) = z in ^7, 

(3.46) { div(eVx* - e-f{w*j) = in ft, 

ty* = 0, X* = on r, 
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which is the unique solution of 

(3.47) a*(K,x*),K,x'))- I z-w'dx,\l{w',x')^V, 
where 

a*{{w, x), {w',x')) = I {{<y{w) + e^Vx) : l{w') + (eVx - e-i{w)) ■ Vx'} dx, V(w', x') e V. 

Again this problem has a unique solution since the bilinear form a* is also coercive on V . 
Since the system (I3.46P is strongly elliptic, we deduce that (w*,x*) belongs to H^[VlY' x 
H'^{n) with the estimate (see Theorem 10.5 of [IJ or Theorem 4.5.3 of [H|) 

(3.48) \\w*\\2,n + \\x*h,n<C\\z\\o.n, 

where here and below C is a positive constant that depends only on aijki,£, lJ-,eijk and on 

n. 

By using the differential equations from (j3.46p . we may write 

z\'^dx = j V{(T{w*) + e^Vx*) ■ zdx 
n Jn 

{S/{a{w*) + e Vx*) • z + div(£Vx* — e7(w*))x} dx. 
n 

Applying Green's formula we get 

z\^ dx = - {{(^(w*) + e^Vx*) : 7(z) + (eVx* - ej{w*)) ■ Vx} dx 
Q Jn 

+ I {(j{w*) + e^Vx*)i'-zdS 

' {{a{z) - e^Vx) : l{w*) + (eVx + 67(2)) • Vx*}dx 
n 



+ {a{w*) + e'Vx*>-zdS. 

Applying again Green's formula and reminding problem p.40[) . we have found that 

|zpda;= f{a{w*)+e'^Vx*>-udS. 



By Cauchy-Schwarz's inequality and the estimate (|3.48p (with the help of a trace theorem), 
we obtain finally 



z\^dx<C / \u\'dS. 
n Jr 



This proves f^AT\ because 4 /p |up dS* < £(i). 

By deriving the system p.40p in time, the estimate (|3.4ip also shows that 



/ \dtz\^dx<C f \dtu\^dS. 
Jn Jt 



This yields (|02ll owing to the identity (|2?25l) . ■ 

At this stage we need to exploit the fact that eVx + 67(2) is divergence free, hence it is 
the curl of -0 S Xt{^), where 

Xt{^) = {(/)£ H^{Q.f : divV' = in f7, and ^' • i^ = on T}. 

More precisely we have the following result. 
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Lemma 3.6 Let {u,E,H) be a strong solution of U.7\) and {z,x) G H^{fl)'^ x Hq{Q) the 
weak solution of ^3.40^ - Then there exists ip G Xt{^) such that 

(3.49) eVx + 67(2) = curl?/;, 
with the estimates 

Of 

(3.50) \\^\\l^<c\\u\\l^<—S{t), 

(3.51) l|5tV'llo,o < C\\dtu\\l^ < -Cdt£{t). 
where C is a positive constant independent of t. 

Proof: We remark that (see (|3.40p ) eVx + 67(2) is divergence free in fi, hence as Vt is 
simply connected, we deduce (see Theorem I. 3.5 in [llj) that there exists ■0 G Xt{^) such 
that (|3.49p holds with the estimate 

||0||i,o<C||eVx + e7(z)||i,o. 
Thanks to (1X15)) . we get 

(3.52) \mi.n<C\\u\\i,n<CE{tY/\ 

Let us finally consider the problem: find x solution of 

{curlcurlx = f/; in $7, 
divx = in 17, 
X- V — {)^ curlx X z^ — on F. 

The variational formulation of this problem is: find x G ^t(^) solution of 

(3.54) h{x,0) = I i:-edx,ye e Hrin), 

Jn 

Kx, S)^ I {curlxcurl6' + divxdiv6'} dx, Vx, 9 G Ht{^), 



where 



and 



Hrin) = {<p£ H\nf : (^ • 1/ = on F}. 



It is well known (see for instance [7J) that b is coercive on Ht{^) and therefore problem 
(|3.54p is well posed, its solution x furthermore satisfies (|3.53p because V' is divergence free. 
Moreover as the system curlcurl— Vdiv — —A is strongly elliptic and the boundary conditions 
in (|3.53p cover this system, we get that x belongs to H^{il)^ with (see again Theorem 10.5 
of [IJ or Theorem 4.5.3 of \^) 

(3.55) llxll2,o<C||V^||o,n. 

Now as before we can write by using Green's formula and the identity (|3.49p 



2 



tp ■ curlcurlx dx 
n 

curlf/; . curlx dx 
n 



/ (eVx + 67(2)) • curlx (ia; 
Jn 



/ V(e curlx) • z da; + / (e curlx)i^ • ^^rfS'. 
Jn Jr 
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By the estimate p.55p and reminding that z = u on F, we obtain 

||^||o,n<C(||z||o,n + |k||o,r). 
By the estimate p.4ip . we arrive at 

\mln<C\\u\\lr, 
and we conclude as in the previous Lemma. ■ 

Lemma 3.7 Let {u,E,H) be a strong solution of 111.7]) . Then for all 9 > there exists 
a constant C{9) > (which does not depend on T hut depends on 9, the domain and the 
coefficients aijki,e, fi,eijk, A) such that 

(3.56) / \u\^dSdt<C{9)£{0)+9 f £{t)dt. 



Proof: We multiply the first identity of (|1.7p by z e H^{Vl)^ from Lemma 1531 and integrate 
on Qt to get 

z ■ {dtu - \7a{u, E)) dxdt = 0. 



By Green's formula we obtain 

(z • d'^u + (t(w, E) : 7(z)) dxdt - / z ■ {cr{u, E) ■ v) dSdt = 0. 

Qt .JJ2t 

Using the second boundary condition in (ll.7p and the boundary condition in p.40p . we 
obtain 

A I \u\^dSdt^~f u-{dtU + Q{E xi'))dSdt- I {z ■ d^u + a(u,E) : -f{z))dxdt. 



Owing to (|3.44p we arrive at 



A I \u\^dSdt< 



u ■ [dtu + Q{E X v)) dSdt - / (z • dfu + e^Vx : 7(m) - e^S : j{z)) dxdt. 

St ^ Qt 

By using the identity ej{u) = D — eE, we get 

(3.57) Af \u\^dSdt<~ I u-{dtu + Q[Exv))dSdt 

- / (z-dtU + Vx-D-E- (e7(z) + eVx)) dxdt. 
Jqt ^ ' 

We now transform the two last terms of this identity, first by Green's formula in space, we 
see that 



Vx • E dxdt = xdivZ? dxdt + / xD ■ v dSdt = 0, 

J Qt "^ 5jt 

since D is divergence free and x = on F. On the other hand, by the identity (|3.49p we 
have 

/ E ■ {e"f[z) + eVx)dxdt = / E ■ cuyI^j dxdt, 
Jqt Jqt 
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and by Green's formula in space 



E ■ (e7(z) + eVx) dxdt = / curl£; • ip dxdt + {E x ly) ■ -ijj dSdt. 
Now reminding that iidtH = curl_B and using an integration by parts in time, we arrive at 

/ E ■ {ej{z)+e\7x)dxdt^ fiH ■ ^t^p dxdt + fj.H ■ Tp dx + {E x u) ■ ipdSdt. 

Jqt Jqt Jn Jst 

In the same manner an integration by parts in time yields 



z ■ dfU dxdt = I dtz ■ dtu dxdt + j z ■ dtu dx 



These identities in p.57p lead to 



(3.58) A I \u\^dSdt< - I 



{u ■ {dtu + Q{E X v)) + iE xiy) -tp) dSdt 



+ j {dtzdtu + ^iH ■ dtip) dxdt — z ■ dtu dx 

Jn 



^iH ■ ip dx 



It remains to estimate each term of this right-hand side. For the first term applying 
successively Cauchy-Schwarz's inequality, Young's inequality and the identity (|2.25p we may 
write 



u-{dtu + Q{Exv))dSdt 



: 4/ \u?dSdt+^ f {\dtu\'^ + \Exv\^) dSdt 
Since the energy is non-negative, we arrive at 



(3.59) 



u ■ dtu dSdt 



<^J^ \u\'dSdt + ^£{0). 



For the second term by using Cauchy-Schwarz's inequality, Young's inequality, a trace 
theorem, the estimate p.52p and again the identity (|2.25p 



{E xiy)-ij dSdt 



< 



< 



< 



d f \mlndt+^ [ \Ex,y\^dSdt 

Jo '' JEt 

e f £{t)dt + ^ f \Ex v\'^dSdt 
Jo c* Jj^^ 

f £{t) dt-^ f dt£{t) dt. 
Jo ^ Jo 



As before the energy being non-negative, we arrive at 



(3.60) 



{Exv)-ij; dSdt 



<e I £{t)dt+^£iO). 
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For the third term we use successively Cauchy-Schwarz's inequality, Young's inequality, 
the estimate p.42p and the definition of the energy to get for all 6* > 



dfZ ■ dfudxdt 



< — / \dtz\^ dxdt + - I \dtu\^dxdt 
2c^ Jqt 2 Jq^ 



< 



C_ 
20 



dt£{t) dt + e £{t) dt. 
Jo 



Again we get 

(3.61) / dfZ ■ dtudxdt 

Jqt 

As for the third term replacing the estimate (|3.42l) by p.5ip we get for the fourth term 



< ^£{Q) + e I £{t)dt. 





(3.62) 



^.H ■ dt^ dxdt 



<^£i0) + 9 f £{t)dt. 
9 Jo 



For the fifth term the application of Cauchy-Schwarz's inequality, the estimate p.4ip and 
the definition of the energy directly gives 



(3.63) 



z ■ dtudx 



< C{£{0) + £{T)) <2C£{0) 



since the energy is non-decreasing. 

Similarly using (|3.50p instead of p.4ip . we have 



(3.64) 



fiH ■ il) dx 



< Cf (0). 



The estimates (|3.59p to p.64p into the estimate p.58p yield the conclusion. 
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